Preprint typeset in JHEP style - PAPER VERSION 



HPro: A NLO Monte-Carlo for Higgs production via 
gluon fusion with finite heavy quark masses 



Charalampos Anastasiou 

Institute for Theoretical Physics, ETH Zurich, 
8093 Zurich, Switzerland 
E-mail: babis@phys.ethz.ch 

Stefan Bucherer 

Institute for Theoretical Physics, ETH Zurich, 

8093 Zurich, Switzerland 

E-mail: stef abuOitp .phys . ethz . ch 

Zoltan Kunszt 

Institute for Theoretical Physics, ETH Zurich, 

8093 Zurich, Switzerland 

E-mail: kunsztOitp . phys . ethz . ch 



Abstract: We compute fully differential next-to-leading order QCD cross-sections for 
Higgs boson production via gluon fusion in the Standard Model. We maintain the full 
dependence of the cross-sections on the top and bottom quark mass. We find that finite 
quark mass effects are important given the achieved precision of QCD predictions for 
gluon fusion. Our Monte-Carlo program HPro can correct existing NNLO fully differential 
calculations, which employ the approximation of an infinitely heavy top and a vanishing 
bottom quark Yukawa coupling, for heavy quark finite mass effects through NLO. 
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1. Introduction 



QCD radiative effects in the gluon fusion process for the production of Higgs bosons at 
hadron colliders have been the subject of detailed theoretical investigations in the last two 
decades. The next-to-leading order QCD corrections to the inclusive cross-section were 
computed already in the nineties [1-3]. In these works it was shown that QCD perturbative 
corrections are substantial. 

Exact calculations of perturbative corrections in the gluon fusion process, such as in 
Ref. [3] , are technically involved due to the presence of massive quark loops already at the 
leading order. With these pioneering NLO computations, the quality of the simplifying 
approximation of an infinitely heavy top quark and vanishing Yukawa couplings for all 
other quarks could be assessed. 

Such an approximation is indispensable for computing QCD perturbative corrections 
beyond NLO. Next-to-next-to- leading order QCD corrections for the gg — > H inclusive 
cross-section were computed in Refs [4-6]. NNLO corrections, within the same approxima- 
tion, for fully differential cross-sections were computed in Refs [7-9] and in Refs [10,11]. 

The theoretical uncertainty due to scale variations of NNLO inclusive and differential 
cross-sections is by now remarkably small. While existing tools, such as HIGLU [12] can 
be used to correct the predictions for the total cross-section for finite quark-mass effects 
through NLO, fully differential cross-sections for the Higgs boson and its decay products 
cannot be corrected for the same effects with existing tools. We remedy this situation in 
our publication. 

We have written a parton level Monte-Carlo program HPro which computes fully dif- 
ferential cross-sections at NLO in QCD while keeping the exact dependence on the finite 
top and bottom quark mass. HPro includes the decays of the Higgs boson to photons and 
four-lepton final states. It can be used in conjunction with the fully differential NNLO pro- 
gram FEHiP [8,9], correcting for finite quark mass effects through NLO. An earlier version 
of HPro has been used in [13] to estimate the finite bottom mass effects. Later in a similar 
study [14] these effects have been also accounted for using HIGLU. 

We present the calculation method in Section 2. In Section 3 we use HPro to compute 
the NLO total cross-section and review the quality of the approximations which are usually 
made for heavy quark loops in the gluon fusion process. In Section 4, we compute for 
illustration various kinematic distributions and study the effect of finite quark masses to 
their shape. We present our conclusions in Section 5. 

2. Implementation 

The computation of the NLO corrections to the gluon fusion process requires the two- 
loop gg — > h amplitude. This has been first computed by means of an one-dimensional 
integral representation derived in Ref. [3]. In Ref. [15], this result was expressed analytically 
in terms of harmonic polylogarithms using a "series expansion and matching" method. 
Independent analytic evaluations of the two-loop amplitude were performed in Ref. [16] 
and in Ref. [17]. 
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At NLO also real radiation sub-processes gg — > gh, qg — ► qh and qq — > g/i contribute. 
The corresponding matrix elements have been first computed in Ref. [18]. In Ref. [19] 
these matrix elements have been expressed in terms of standard one-loop scalar integrals 
and helicity amplitudes. We have recomputed these contributions along the lines of Ref. [19] 
and found full agreement. In particular we have also compared to a very recent calculation 
in Ref. [20]. For completeness we present the results in Appendix A. 

Virtual and real corrections develop singularities which only cancel in their combination 
and by adding the appropriate collinear counter-term for PDF evolution. In order to 
obtain an expression suited for numerical integration we apply the FKS subtraction method 
[21]. Recently this method has been applied in Ref. [22] for the matching of NLO Higgs 
production to shower Monte Carlo in the heavy top mass approximation. 

We note that in the case of total cross sections the straightforward application of the 
FKS subtraction method leads to the same analytic formula for the subtraction terms 
that has been found in Ref. [3]. With appropriate insertion of measurement functions we 
easily get analytic expressions suited for the calculation of any differential distribution. In 
practice one might think of the measurement function as a vector valued function with 
each component being a bin of a certain distribution. There exist implementations of the 
VEGAS algorithm supporting vector functions, such as the VEGAS routine in the CUBA 
library [23], which allows to obtain a reliable Monte-Carlo error estimate for each bin. In 
combination with dedicated decay routines for the diphoton and leptonic final states, this 
feature enables us to produce various distributions relevant for experimental searches of 
the Higgs boson, corrected by finite mass effects, as we will demonstrate in the following 
sections. 

3. Total cross-section at NLO 

In this section, we revisit the gluon fusion cross-section at NLO. This serves as a check of our 
Monte-Carlo HPro against the predictions of HIGLU [12], and to emphasize the importance 
of finite quark mass effects in Higgs boson production. For the numerical results of this 
paper we use MSTW 2008 parton distribution functions [24]. 

We begin our study by revisiting the total cross-section in the LO and NLO approxi- 
mation as a function of the Higgs boson mass (Fig. 1). As it is well known [1-3], NLO QCD 
perturbative corrections are substantial. We note here that the perturbative corrections 
are slightly smaller with the latest parton densities [24], mainly due to the higher value 
of a s used at leading order. We also note that a significant QCD correction is found at 
NNLO [4-6], which is not included in the results of this article. NNLO corrections stabilize 
the perturbative expansion and reduce the scale variation to the ~ 10% level. However, 
NNLO computations rely on an approximate treatment of heavy quark loops. 

In Fig. 1, we show the effects of different treatments for the heavy quarks. With HPro, 
we compute the exact LO and NLO cross-sections, where all loop diagrams with massive top 
and bottom quarks are evaluated exactly (we denote with "top+bottom" the corresponding 
cross-sections in the plots of this paper). An approximation which can be made, is to 
consider a vanishing bottom Yukawa coupling and to evaluate exactly only the top-quark 
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Figure 1: Total Cross-Section at Tevatron and LHC. 



loops; in our plots, we denote this approximation as "top-only". NNLO computations 
are performed in what is known as the infinitely heavy top-quark approximation. In the 
"witop = oo" approximation, the bottom Yukawa coupling is set to zero, and the cross- 
section at higher orders is estimated by the formula, 

(^top-only ' 
(N)NLO i 
^gpiy-l M> 

where bottom quark loops are ignored, and the leading order cross-section is reweighted 
with the ratio of the cross-sections at higher orders and the leading order in the limit of 
an infinite top-quark mass. In order to compare the different approximations we introduce 

SX i = ^— (3.2) 

where X is the inclusive cross section or a normalized differential cross section and i is 
labelling the contribution (top+bottom, top-only, etc.) 




m h [GeV] m h [GeV] 



Figure 2: Percent differences of the exact NLO (LO) total cross-section with finite top bottom 
masses or the NLO (LO) total cross-section with exact top mass effects but zero bottom Yukawa 
coupling with respect to the usual approximation at Tevatron and the LHC. 
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In Fig. 2, we show the cross-section deviations from the m top = oo approximation of 
Eq. 3.1 when "top-only" (blue) and the complete "top-bottom" mass effects are taken into 
account. "Top-only" contributions are approximated within a couple of a percent up to 
the mn = 2mt op threshold. However, for a light Higgs boson, bottom quark contributions 
are important and can reach ~ —8%. It is then important that bottom loops are taken 
into account for a precise evaluation of the total cross-section [13]. We also observe that 
the contribution from bottom-quark loops decreases at NLO in comparison to LO. 
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Figure 3: Ratio in percentage of top-only, bottom-only and top-bottom interference components 
with respect to the total cross-section at NLO, for Tevatron and LHC. 

In Fig. 3, we show the relative contributions to the NLO total cross-section from top- 
quark loops only, bottom quark-loops only, and from the interference of top and bottom 
loops. Top-only contributions are dominant, while bottom-only contributions are negligible 
over the whole Higgs mass range. Top-bottom interference terms are important at the few 
percent level and are negative for a light Higgs boson. It should be noted that the relative 
importance of the three contributions for a heavy Higgs boson or a pseudo-scalar Higgs 
boson in the MSSM may be drastically different than in the Standard Model [25]. 

The magnitude of QCD corrections depends strongly on the mass value of the heavy 
quark in the loops if we use the pole mass for the qqH coupling. In Fig. 4 we consider the 
cross-section for the gluon fusion cross-section at the Tevatron and the LHC, considering 
only one heavy quark with a mass m g . We find that, in the pole scheme, the K-factor 
is reduced significantly for small values of the quark mass consistent with neglecting the 
running of the quark mass in the Higgs coupling at leading order. 

For the Standard Model where both top and bottom quark loops contribute to the 
gluon fusion process, we study separately the magnitude of QCD corrections for the top- 
only, bottom-only, and top-bottom interference terms. The corresponding K-factors at 
NLO are plotted in Fig. 5. While the top-only contributions receive a large K-factor, NLO 
QCD corrections to the top-bottom interference and bottom-only terms are milder. As a 
consequence, the importance of the bottom-quark loops is smaller at NLO than at LO. 

The small NLO QCD corrections to the top-bottom interference contribution, which 
is also a very small fraction of the top-only contribution, suggests that a more precise eval- 
uation at NNLO is not necessary. The top-only contribution receives however large NLO 
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Figure 4: NLO cross-section and K-factor for the gluon fusion process via one only heavy quark 
at the Tevatron and the LHC, as a function of the quark mass. 
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Figure 5: NLO K-factors for the top-only, bottom-only, and top-bottom interference contributions. 
We restrict to a range of toh where top-bottom interference contributions are still sizable. 



corrections and it requires an evaluation at NNLO. As shown in Fig. 2, this contribution 
can be approximated using Eq. 3.1 better than 2% for a light Higgs boson, and better than 
10% for a Higgs boson with a mass above the top-pair threshold. It appears to us, that the 
combination of the NLO cross-section with full dependence on the top and bottom quark 
masses and the NNLO correction using the approximation of Eq. 3.1 yields a very precise 
estimate of the gluon fusion cross-section, where differences with an NNLO calculation 
with exact finite quark mass effects should be quite small. 
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The results of this section have been extensively cross-checked with HIGLU and excellent 
agreement has been found. 

4. Differential cross-sections at NLO with finite quark masses 

The search for a Higgs boson at hadron colliders is complicated due to the large cross- 
sections of background processes. Sophisticated experimental analyses are required, where 
it is essential to find optimized selection cuts. In addition, it is often necessary to perform 
a detailed probabilistic comparison of measured shapes for kinematic distributions with 
theoretical predictions for the signal and background processes. The role of very accurate 
Monte-Carlo programs which are fully differential is very important for these purposes. 

The fully differential NNLO Monte-Carlo's, FEHiP [7-9] and HNNLO [10,11], are avail- 
able for the gluon fusion Higgs boson production process. Given the complexity of NNLO 
computations, these programs employ the approximation of Eq. 3.1. In some cases, ex- 
perimental cuts lead to a significantly smaller scale variation than in the total cross- 
section. This enhances the importance of other uncertainties, such as the one due to un- 
accounted finite quark mass effects. A characteristic example is the accepted cross-section 
for pp — > H — > WW — > llvu where a jet-veto and other cuts reduce the uncertainty due to 
scale variations by a factor of about two with respect to the total cross-section [9, 11,28]. 

With our exact NLO Monte-Carlo HPro, we can correct the predictions of FEHiP and 
HNNLO for finite quark mass effects through NLO. HPro can be used to compute fully 
differential cross-sections and distributions at NLO for the Higgs boson and the final state 
particles in the two photon and four lepton decays. HPro is merged with the FEHiP NNLO 
Monte-Carlo and it will be released in a forthcoming publication. In this section, we 
illustrate the shapes of a few kinematic distributions with our new HPro NLO Monte-Carlo 
and compare them with the corresponding predictions in the "mtop = oo" approximation. 
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Figure 6: NLO normalized rapidity distribution at the LHC, tor = 150 GeV. Finite quark mass 
effects affect the shape of the distributions only in the very high rapidity range, where only a tiny 
fraction of events take place. For y H values larger than the one marked by the grey line, less than 
10 -3 of the total number of events take place. 

In Figs. 6 and 7, we study the normalized NLO rapidity distribution of a Higgs boson 
with mass mn = 150 GeV at the LHC and the Tevatron, respectively. Except for very large 
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Figure 7: NLO normalized rapidity distribution at the Tevatron, mn = 150 GeV. Finite quark 
mass effects affect the shape of the distributions only in the very high rapidity range, where only a 
tiny fraction of events take place. For y H values larger than the one marked by the grey line, less 
than 10 -3 of the total number of events take place. 

rapidities, where almost no events occur, the distribution is not affected by mass effects. 
Even in that range the deviations from the "mtop = oo" limit are less than 5% at the LHC 
and less than 10% at the Tevatron. We conclude that small- x effects are moderate and 
without phenomenological consequences in the case of the rapidity distribution at NLO 
shedding light on an open question raised in [26,27]. 

Another important differential distribution is the transverse momentum of the Higgs 
boson. Finite quark mass effects for the pr distribution have also been studied in earlier 
publications [18,19,25], and recently both electroweak and finite quark mass corrections 
were computed and combined [20]. 




Figure 8: Normalized px distribution at Tevatron, toh = 120 GeV. Massive corrections are 
important for large pt- However only a very small fraction of events exists in this range. 

In Fig. 8 (left panel) we present the normalized cross-section at the Tevatron in pi 
bins of 2 GeV, for a Higgs boson with mass mn = 120 GeV. At small values of pi the 
bin cross-sections cannot be computed accurately in perturbation theory, see Fig. 9, and 
an all orders resummation is required [29,30]. A meaningful result is obtained, however, 
when the bins at low pi are added up together. In order to study the effect of finite quark 
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Figure 9: First few bins of pt distribution at Tevatron, toh = 120 GeV. The values of these bins 
are unphysical and require resummation. 

masses it is more convenient if we demonstrate uncombined low pt bins. For this purpose, 
we present the pt distribution in the approximation of Eq. 3.1 (cyan), in the "top-only" 
approximation (blue) where the bottom loops are ignored but the top-loops are evaluated 
exactly, and with the complete "top-bottom" mass dependence (red). We have compared 
our results with the authors of [20] and found full agreement within numerical errors. In 
the right panel of Fig. 8 we show the percent deviations of the complete result (red) and the 
"top-only" approximation for the normalized p<r distribution from the approximation of 
Eq. 3.1. At small pr, there are very small differences due to finite quark-mass effects. We 
observe some important shape deviations due to the effect of top and bottom quark loops 
at intermediate pr- As it has already been observed in Ref. [20] finite quark effects are 
very large at high pi, where the quark production channel becomes dominant. Note that 
additional electroweak corrections affect the shape considerably [20]. At a higher Higgs 




Figure 10: Normalized pt distribution at Tevatron, ttih = 180 GeV. Compared to hih = 120 GeV, 
mass corrections lead to smaller deviations from the "mtop = oo" approximation. The spectrum of 
the bottom-only contribution is much softer (green, left panel). 

boson mass value of mn = 180 GeV (Fig. 10) we find an even milder effect at low pt, 
while the magnitude of the deviations at a large px is somewhat reduced but still large. 
For phenomenological purposes, these large deviations concern a tiny fraction of potential 
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Higgs signal events for both me = 120 GeV and mn = 180 GeV mass values. We show in 
Fig. 10 also the "bottom-only" contribution and observe that in this case the pt spectrum 
is much softer, as pointed out in Ref. [25]. 




Figure 11: Normalized px distribution at LHC for win = 120 GeV and toh = 180 GeV. Mass 
corrections are much more modest than at Tevatron where the qq channel plays a much bigger role. 

At LHC energy (14TeV), Fig. 11, we observe significant bottom-loop effects for a light 
Higgs boson (mn = 120 GeV). These are reduced, for a heavier Higgs boson with mass 
run = 180 GeV. Shape deviations due to finite quark mass effects can reach up to 10% 
at high px- It is interesting that bottom quark loops for a light Higgs boson change the 
shape at low px- As we explained, the fixed order px spectrum is not physical at low px- 
However, these deviations may also survive after a complete resummation is performed via 
the matching procedure. It is interesting to examine the px spectrum for a Higgs mass 
where the heavy top approximation is formally invalid. In Fig. 12 we plot the normalized 
distribution for mn = 400 GeV at the LHC. Deviations of the "top-only" contributions 
from the infinitely heavy top-quark approximation are small for px < 80 GeV. At higher 
Pt the difference increases. 

Finally, we present normalized distributions for Higgs decay final state. In Fig. 13 we 
present the pseudorapidity difference and average pt distribution of the two photons in 
the process pp — > H — > 77. Finite quark-mass effects do not affect these distributions. At 
higher Higgs boson masses the process pp — * H — > WW — > llvv is dominating and we show 
as an example for this decay mode the 4>u distribution in Fig. 14. <f>n is the angle in the 
transverse plane between the two charged leptons in the final state and we find again, that 
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Figure 12: Normalized pt distribtution at LHC for to,h = 400 GeV. Mass effects become important 
as TOtop = oo approximation is formally invalid. 
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Figure 13: On the left: Normalized distribution of the average photon transverse momentum, 
p T avg _ (p T 7i -fp T 72 j j2 w On the right: Normalized distribution of photon pseudorapidity difference, 
y* = \rpi — 7y 72 |/2. In both plots mn = 120 GeV and we assume LHC energies. 
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the shape is very well reproduced by the "mt op = oo" approximation. 

In summary, we have found that the shapes of distributions for leptons and photons 
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from the decay of a Higgs boson are very well approximated by Eq. 3.1. In addition, 
accepted cross-sections after the application of cuts on jets are affected consistently with 
the expectations from the shapes of the Higgs pt spectra. 



5. Conclusions 

In this paper, we have presented a NLO partonic level Monte-Carlo program HPro, which 
computes the top and bottom quark mass dependence of differential cross-sections exactly. 
HPro computes accepted cross-sections after selection cuts and kinematic distributions for 
final-state particles in the diphoton and four-lepton decay channels. 

Finite quark mass effects are important and can affect the precision of NNLO calcu- 
lations in the infinite top-quark mass approximation. We can use HPro to correct NNLO 
differential cross-sections for finite quark mass effects. For this purpose, we have interfaced 
HPro with the fully exclusive NNLO Monte-Carlo FEHiP. The combined program will be 
released in a forthcoming publication. 
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A. Real radiation matrix elements 

Here we present the results for the matrix elements squared for real radiation processes 
gg — > gh, qg — > qh and qq — > gh in terms of helicity amplitudes. We introduce the following 
shorthand notation for scalar one-loop integrals, 



where Bq, Cq and -Do are the standard one-loop integrals in the notation of e.g. [31]. B { 
is the finite part of the Bq function, i.e. 



There exist several publicly available packages for evaluating one-loop integrals [31-33]. 
However, we used a private implementation and checked against these packages. 

In the following we present the real radiation matrix elements squared, averaged over 
spin and color and divided by the flux factor, in terms of helicity amplitudes, closely 
adapting the notation in [19]. Note however the different convention regarding helicity 



B (s, m) 
C (s, m) 
Ci (s, t, m) 
D (s, t, m) 




(A.l) 
(A.2) 
(A.3) 
(A.4) 




(A.5) 
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labels. The simplest process is qq — > 5/1 and the result can be expressed by a single 
independent helicity amplitude, 



T^-( S12 , Sl3 , * 23 ) = "[^N*^ ( I £ Mq(s\2, ,13, *23) Mq(s\2, *23, -is) 



with 



(A.6) 



^Q 9 (S12,S13,S 23 ) 



2(si 3 + s 23 - 4mg)Ci (m^, s i2 , m Q ) 



Aq"IqS23 
S 2 3 + S13 

H ^ — (B (s 12 ,m Q ) - B (m#, m Q )) - 4 

«23 + s 13 



(A.7) 



Here, tuqAq (Aq = 1/v), is the Higgs-quark-quark coupling and we sum over heavy quarks, 
Q = t,b. 

Similarly simple is the result for qg — > qh and — > qh sub-processes which is obtained 
by crossing, 



Tg ff (si 2 ,si 3 ,s 23 ) 
T^(si 2 ,si 3 ,s 23 ) 



AT, 



N 2 _ 1 -qq 



T^(s 23 ,si 3 ,si 2 ) , 



7V C 2 -1 qq 



YL(S13,S1 2 ,S 23 ) 



(A.8) 
(A.9) 



More involved is the expression for the gg — > <?/t process. However, it has a compact 
representation in terms of only two independent helicity amplitudes, 



Tp„(si 2 ,Si 3 ,S 23 ) = 



a s 3 N c 



E M q!++-( S 12>S13,S23) 



87r(iV c 2 - l)s 2 12 s 13 s 2 3 

+ | E M q! ++ -( S 13, 512, ^23)| 2 +| E M Q; ++ J^3, *13, Surf + \ £ M 9 Q 9 . +++ (s 12 , S 13 , 823^ ■ 

(A.10) 
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The amplitudes appearing in this expression are given by 

4si2(>i2 2 - S13S23) 



Mq? ++ _(si2, S13, S23) = Agm^ 

1 Sl2S23(4m|si 3 - S12S13) 

" 2 ^ 

1 S12S13(4S23"T-Q - S12S23) 
" 2 



(«23 + Sl 2 )(si3 + S12) 
D(S12,S23,?^Q) 

D (si 3 ,si2,mQ) 



1 S23Sl3(-Sl2 2 + I2s 12 m% + 4si 3 S23) 

+ 77 D(S13,S23,"IQJ 

^ «12 

4si 3 (2si 2 S23 + «23 2 ) , \ ^S 2 3{2s 12 S 13 + S 13 2 ) 

- f : T2 B ( s i3, nKj 7 7 B ( s 23, mq) 

(S23 + Sl 2 )^ (Sl 3 + Si 2 )^ 

- 7 ; 2 to 3 / 23 , ( ~ 4si 3 2 Si2 - 2s 2 3Sl3 2 ~ 8si 3 Si 2 S 2 3 

(Sl3 + Sl2j (S23 + Sl2j V 

- IOS13S12 2 - 2s 2 3 2 si3 - IOS23S12 2 - 4si 2 S23 2 - 8si 2 3 )b (m 2 H ,mQ) 

(S13«12S23 - 4si 3 S23m^)(s 2 3 + S13) 

H Ci {m H ,si2,m Q ) 

S13S23 

+ 7 7 r(4si 3 S23 2 - S12S23 2 + 2 S2S + 2si 3 S 2 3Sl2 + «12 3 

(S23 + S12) S12 

+ 4mg(s 2 3 2 + 2si 2 s 2 3 - si2 2 ))Ci (m 2 H , s 13 ,m Q ) 

1 (a 2 2 , S 13 3 S23 . ,3 
(4s 2 3«13 - S12S13 + 2 h 2S23«13S12 + «12 



+ Am 2 Q (s 13 2 + 2si 2 si3 - si2 2 ))Ci (m 2 H , s 23 ,m Q ) 

~Si 3 2 S 23 , \ S23 2 Sl3 r1 / \ 

2 C (si 3 , m Q ) - 2 C (s 2 3, m Q ) 

S\2 S\2 



and 

A*Q? +++ (si2, S13, S 23 ) = AqTOq 



1 523512(4^,^513 


- S13S23 


- Sl2«13 


" *13 2 ) 


2 


S13 






1 S12S13(4S23"IQ 


- S12S23 


- Sl3«23 


" «23 2 ) 


2 


«23 






1 S23Sl3(4si 2 m^ 


- Sl2«23 




S12S13) 



2 ai2 

(S13«12S23 + S13S23 2 + Sl3 2 «23 - 4S1 3 S 2 3"T-^)(S23 + «13 



H - Ci (m H ,s 12 ,mQ) 

■S23S13 

(«12S23 2 + «12 2 S23 + S13S12S23 - 4si2S23m^)(s 2 3 + S12) / 9 \ 

H Ci (m H , si 3 , m Q ) 

■S12S23 

(Sl2 2 «13 - 4si 2 m^Si 3 + S13S12S23 + Sl 2 Sl3 2 )(si3 + S 12 ) , 

H Ci {m H ,s 23 ,m Q ) 

S12S13 



(A.11) 



4(S13 + S23 + Sl2) 

D (S12, S 2 3,"2q) 
D(si3,Si2,mg) 

D(si3,s 2 3,m Q ) (A.12) 
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